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ABSTRACT 

In this paper we discuss the spin-equilibrium of accreting neutron stars in LMXBs. 
We demonstrate that, when combined with a naive spin-up torque, the observed data 
, leads to inferred magnetic fields which are at variance with those of galactic mil- 

■ lisecond radiopulsars. This indicates the need for either additional spin-down torques 
C^ - , (eg. gravitational radiation) or an improved accretion model. We show that a simple 

consistent accretion model can be arrived at by accounting for radiation pressure in 
rapidly accreting systems (above a few percent of the Eddington accretion rate). In 
I ■ our model the inner disk region is thick and significantly sub-Kcplerian, and the es- 

t^J" ! timated equilibrium periods are such that the LMXB neutron stars have properties 

■ that accord well with the galactic millisecond radiopulsar sample. The implications 
^ (— I for future gravitational-wave observations are also discussed briefly. 

8 ■ 1 INTRODUCTION 

£3 ■ 

^ ' In the last few years the evidence in favour of the notion that neutron stars are spun up to millisecond periods in accreting 
• • ■ systems has strengthened significantly. The discovery of the millisecond X-ray pulsar SAX J1808. 4-3658 in an accreting 
, ^ ' low-mass X-ray binary (LMXB) provided the long anticipated missing link between the general LMXB populatio n and the 

■ millisecond radio pulsars. Since then four similar system has been observed, further strengthening the connection iWiinandsl 
?H ' 120031) . Furthermore, the link between the twin-peak separation of the kHz quasiperiodic oscillations (QPOs) seen in a number 

of systems and the spin of the neutron star has become somewhat clearer (although the underlying mechanism is still under 
debate) following the o bservation of QPOs in SAX J1808. 4-3658. It now appears as if the QPO separation is either equal to 
or half the spin period <|Millerll2003h . 

When the first indications of rapidly spinning neutron stars in LMXBs were discussed more than five years ago, the 
results suggested that the systems were clustered in a surprisingly narrow range of spin frequencies 250-370 Hz. As such spin 
rates are far below the predicted break-up limit of about 1 kHz, the data pointed towards the presence of a mechanism that 
could counteract the accretion spin-up torque. The ob vious candidate — the interaction between the accretion disk and the 
magnetosphere of the neutron star — was discussed bv I White fc Zhanel l)l997t ). Their results seemed to indicate the need for 
an unanticipated link be tween the accret ion rate and the magnetic field strength. Since there is no reason to expect such fine- 
tuning in these systems, iBildsterl Jl998l^ argued that an additional spin-down mechanism may be in operation. He proposed 
that this torque could be provided by gravitational-wave emission, and th at the required asymmetries wou ld be induced in 
the neutron star crust by accretion. This idea echoed earlier suggestions bv lPapaloizou fc Pringlel Jl978l) and lWaeonerl Jl984l) 
that neutron stars may reach a spin-equilibrium with gravitational waves balancing the accretion torque. 

The possibility that accreting neutron stars may radiate gravitational waves is of great interest given the generation of 
groundbased interferometers (LIGO, GEO600, TAMA300 and VIRGO) that is now reaching design sensitivity. It has been 
recognized that there are three distinct mechanism s that may be able to generate gravitational waves at the required rate. First 
of all, a more detailed study bv lUshomirskv et alJ J2OO0I) suggests that the accretion induced crustal asymmetry proposed by 
Bildsten remains viable. The second possibility is that the stars sp in fast enough tha t the gravitational- wav e driven instability 
of the r-mode oscillations in the neutron star core is activated (see lAnderssonl l|2003h for references). Finally. Icutlerl l|2002h has 
suggested that an internal toroidal magnetic field could lead to unstable free precession resulting in the star "tipping over" 
and becoming an orthogonal rotator, an efficient gravitational-wave source. 
The present investigation is motivated by the following facts: 
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• The observational data has improved considerably since the original discussions in 1997-98. We now know that the LMXBs 
are not clustered in as narrow a range of spins as was originally thought, the current range being 250-620 Hz. It is relevant to 
ask to what extent the more recent data supports the need for an additional spin-down torque, eg. gravitational radiation, in 
these systems. 

• A question that does not seem to have attracted much interest concerns whether a more refined model of the interaction 
between the accretion disk and the magnetosphere of the neutron star would be able to provide a sa tisfactory description o f 
the LMXBs. After all, many important physical mechanisms were not accounted for in the analysis of lWhite fc Zhana dl997^ 
and it may be wise not to refine the various gravitational-wave scenarios before their relevance is investigated. 

• If we suppose that the LMXBs radiate gravitational waves at a significant level, then we need to address many difficult 
issues associated with the detection of such signals. A key issue concerns the spin-evolution of the system. Can we assume 
that the spin-period remains stable on a time-scale of a few months? After all, the signal needs to be integrated for at least 
two weeks in orde r to be detectable in the noisy data-stream. If the system tends to wander, as the data for slower spinning 
systems suggests t'Bild sten et aljll997f) . then we need to be able to model the accretion torque reliably. 

In this paper we aim to address the second of these points. We discuss the argument that an additional spin-down 
torque is needed in the LMXBs, and provide a more detailed accretion model that is able to describe these systems without 
particular fine-tuning of the magnetic field. From this exercise we conclude that it may not be appropriate to assume that 
the neutron stars in LMXBs radiate gravitational waves at a rate that exactly balances the accretion spin-up torque expected 
for a non-magnetic star. We do not think this should be taken as meaning that these systems are irrelevant for gravitational- 
wave physics. The proposed mechanisms for generating gravitational radiation should certainly still work. Yet, our discussion 
makes it clear that modelling these systems is significantly more difficult than has been assumed so far (at least in the 
gravitational- wave community) . Of course, by constructing a more detailed accretion model, we are beginning to address this 
issue. 



2 LMXBS AND THE "STANDARD" ACCRETION MODEL 

In the simplest models of accreting non-magnetic stars it is assumed that matter falling onto the surface of the star provides 
a torque proportional to the angular momentum associated with a Keplerian orbit at the stars equator; 



N « MVGMR (1) 

where M is the mass, R the radius and M the accretion rate. Despite it being well-known that this torque only provides 
an order-of-magnitude estimate, it has been used in most studies of gravitational waves from LMXBs so far. The line of 
reasoning has been that, if the neutron star is at spin equilibrium, then the radiated gravitational waves provide an equal and 
opposite torque. The strength of the gravitational waves can be inferred from the X-ray luminosity, since (assuming that the 
gravitational potential released by the infalling matter is radiated as X-rays) 

provides a link between the observations and the mass accretion rate. 

Accretion onto a magnetised star is different since the pressure of the infalling gas is counteracted by the magnetic 
pressure. By balancing these two pressures (for spherical infall) one obtains the so-called magnetosphere radius 

inside which with the flow of matter is likely to be dominated by the magn etic field. For a s trongly magnetised star the 
magnetic field is expected to channel the accreting matter onto the polar caps iFrank et alJboO^L while the situation may be 
more complex for a weakly magnetised object. 

An approximation of the maximum accretion rate we should expect follows from balancing the pressure due to spherically 
infalling gas to that of the emerging radiation. This leads to the Eddington limit; 

R \ Mr. 



LB x 10- (*) *° (4) 
V10 km/ yr v ' 



yr 

with associated X-ray luminosity 



From these estimates we see that, for accretion at a fraction e of the Eddington rate, eg. M = eMEdd, the magnetic field must 
be accounted for (in the sense that Rm > R) as long as it is stronger than 



AI \ / M 
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(6) 



Since observations indicate that rapidly rotating neutron stars have magnetic fields of the order of 10 s G, and many transient 
LMXBs accrete with e ~ 0.01, we infer that the magnetic field is likely to play a role in these systems. 

The interaction between a geometrically thin disk and the neutron star magnetosphere is a key ingredient in the standard 
model for accretion. The basic picture is that of a rotating magnetised neutron star surrounded by a magnetically threaded 
accretion disk, see Figure [5] for a schematic illustration. In the magnetosphere, accreting matter follows the magnetic field 
lines and gives up angular momentum on reaching the surface, exerting a spin-up torque. The material torque at the inner 
edge of the disk is usually approximated by 



N = MVGMRm (7) 

It is important to note that this torque can be significantly stronger than the rough estimate for non-magnetic stars. 
Meanwhile, outside the co-rotation radius, 



R c « 17 i 



(t^J (omoJ km (8) 

the field lines rotate faster than the local Keplerian speed, resulting in a negative torque. If Rm > Rc the accretion flow will 
be centrifugally inhibited and matter may be ejected from the system. It is easy to see that this will happen if the spin period 
becomes very short, or the rate of flux of material onto the magnetosphere drops. This is known as the propeller regime. As 
accreting matter is flung away from the star in this phase, the star experiences a spin-down torque. To account for this effect 
we alter the material torque according to 



N = MR&[[SIk(Rm) -fi] = MVGMRn 



(9) 



where Q is the spin frequency of the star and and S7k is the angular velocity of a particle in a Keplerian orbit; 

»«(^(^)"" m 

Even though this expression only accounts for the propeller regime in a phenomenological way, it agrees with the expectation 
that accretion will not spin the star up beyond the point Rm = R c - This leads to the equilibrium period 

p nmf B » \ 6/? ( R V 8/ V M V 5/ 7 & V 3/? nu 
P -* °- 3 ° ( ICFG J UtaJ UlMeJ UkJ mS 

Conversely, given an observed spin period we can (assuming that the system is at equilibrium) deduce the neutron star's 
magnetic field. 

Let us now compare this estimate with the observational data, summarised in Tabled To do this we need both stellar 
spin rates and an estimate of accretion rates. Let us first consider spin rates. For the X-ray pulsars, we will use the measured 
pulsar frequency u psl . For those sources that are not pulsars but have burst oscillations, we will assume that the measured 
burst oscillation frequency !A> urs t is the stellar spin frequency. For the third class of sources, which exhibit neither pulsations 
nor burst oscillations, we will use the separation of the kHz QPOs, A^qpo, as an estimate of the stellar spin. As can be seen 
from those pulsars and burst oscillation sources that also have kHz QPOs, this estimate is not exact, as the spin frequency 
can be as much as double the kHz QPO separation. In addition the kHz QPO separation is variable. For the purposes of this 
paper, however, we will estimate stellar spin as the midpoint of the observed range of A^qpo for all sources that do not show 
burst oscillations or pulsations. This places an upper limit on the inferred magnetic field. The inferred field would be lower if 
the true spin rate were greater than Az/qpo- 

The accretion rate can be estimated from the X-ray luminosity, which can be highly variable. It is clear that the estimated 
equilibrium period is shortest when the accretion rate is highest (alternatively for a given spin rate the inferred magnetic field 
is maximal) . In this paper we will assume that the observed spin rate is the equilibrium period associated with the maximum 
accretion rate for a given source, even for sources that are transient or highly variable 1 . This make sense if one assumes that 
the main contribution to the spin-up torque is associated with the phase when the star accretes at the fastest rate. Hence, 
the accretion rates given in Tableware estimates of maximum accretion rates. In addition we will assume R = 10 km and 
M = 1.4M© for all systems. 

Figure ^compares the model's predictions for LMXB magnetic fields with the inferred magnetic fields for the millisecond 
radio pulsars. The agreement is good for LMXBs accreting at the level of 10 _2 MEdd and below. The model does not, however, 
perform well for systems accreting with M w MEdd- The figure shows that the estimated magnetic fields appear to be too 
large for the systems accreting near the Eddington limit, mainly objects for which the spin rate was inferred from the kHz 

1 See however lLamb fc YvJ l2004fl for a discussion of whether neutron stars do reach equilibrium. 
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Table 1. Data for rapidly rotating neutron stars (with spins ab ove 100 Hz), with references given in square brack ets. Source type classifi- 
cations are P (pulsar), A (Atoll), Z (Z source) or U (Unknown) JjTasinger & van der Klisll989tlvan der Klisl20o3) . (T) indicates that the 
source is transient. The frequencies given are pulsar spin frequency (^ psr ), burst oscillation frequency (^b urs t) and separation between the 
two kHz Quasi- Periodic Oscilla tions (Ai/qpo). The accret ion rates shown are estimates of maximum accretion rate , as discussed in the 
main text. References: [11 IWiinands fe van der Klisl <199Sl). [21 [Chakrabartv et alj <2003l). [3] IWiinands et al.l i2003tl . [4] iGallowav et alJ 
]2004l) . [51 iMarkwardt et al.l {20021) [61 iRemillard et al] J2002I). [7 1 IMarkwardt et alJ 120031). [8] C.B.Markwardt , priv ate communica- 
tion, [ 9] IMarkwardt fc Swankl <200a) . [101 IStrohmaver et al] ]2003l) . [Ill IMarkwardt. Swank fc Strohmaverl <2004l). [121 IGallowav et al I 
i2005l) [131 lHartman et alJ i2003l).fl4l iMendez et al.Ul998». [151 D . K.Gal loway, private commun i cation , [161 IKaaret et al] <2002l) [171 
[ Natalucci et alj il99STl . [181 iGiles et alJ <2002l). [ 1 9lTjonker et al.l <2002l). [201 IWiinands et al] l^OOlh. [2 1 ] IZhang et al] <1998h.[22l 
Ismith et al] Jl997ft. ?23l IWiinands fc van der Klisl dl997^ . [241 IKaaret et al] <2003l). [251 IStrohmaver et alJ jl996t), [261 iMigliari et al] 
i2003>). [271 IMarkwardt ,et alj Il999t). [28llFord et al.lT200Cl). [291 IGallowav et alj^OOll). [30llBoirin et a,l\l2000t) . [3!llSmale et al] il98Sl) . 
[32llJonker et _ alj <200d) . [331 IWiinands et al] il99Sl). [3 41 iFord et al.Hl998l) . [351 Ivan Straaten et al] 20031) . [361 1 Jonker et al] <2002lU37 1 
IZhang et al] <1998l). [38 ] IMendez fc van der Klisl feOOOl) . [391 iHoman et al] ]2002l) . [401 iHoman et all ]l999t) . [411 iTomsick et al] Il99sh . 
[421 IZhang et al] <1998l) 



QPO separation. This is, essentially, the conclusion drawn by iBildstenl Jl998l) . There seems to be a need for an additional 
spin-down torque in the systems that accrete at near-Eddington rates. Note that, if we had assumed that the true spin-rate 
was twice A^qpo, as indicated in some of the systems exhibiting burst oscillations, then the inferred magnetic field would be 
roughly half those indicated, which would still be problematic. 



3 A MAGNETICALLY THREADED DISK 

The interaction between an accretion disk and a spinning comp act object involves much poorly known physics. Th e key issues 
were discussed in a number of seminal papers in the late 1970s Jchosh et al l lll977l) : lGhosh fc Lambl Jl978l Il979allbh . see also 
iFrank et alJ J2002I) for an excellent introduction). Although much effort has been invested in this area of research since then 
— after all, accretion is a cornerstone of astrophysics — these early papers remain the "standard" description of the problem. 

In this Section we will focus o n the contribution to t he accretion tor que from a magnetically threa ded, thin disk. Our 
descri ption is based on the work bv lWanel ||1987LIi99, i t|) and lYi et all (jl997h (see also lYi fc Wheeler! Jl99Sl) and lYi fc Grindlavl 
l|l99Sl) ). 

We begin by pointing out that our previous description of the accretion problem was somewhat inconsistent since our 
various estimates, eg., of the size of the magnetosphere, were based on spherical infall of matter. The model can be improved, 
albeit at the cost of introducing several largely unknown parameters. First of all, we need a description of the viscosity in 
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Figure 1. Comparing the neutron stars in LMXBs to the millisecond radio pulsar population. We include all millisecond radiopulsars 
(with periods below 10 ms) in the galaxy. Millisecond pulsar in globular clusters are excl uded since a significant sample of them are 
seen to spin up, an effect likely due to motion relative to the core of the globular cluster .Phinnevl[T993) . which makes the magnetic 
fields inferred for them dubious. In the left panel we compare the inferred magnetic field for the galactic millisecond pulsars, to those 
inferred for accreting neutron stars using the simplest estimate for the spin-equilibrium [B(P cq ) is inferred from Eq. II 11 1. The radio 
pulsars are shown as filled circles, systems showing burst oscillations are represented by open circles, data from systems where the 
spin period is estimated from the kHz QPO separation are open squares and the accreting X-ray pulsars are shown as open triangles. 
We also indicate the (rough) range of magnetic fields for the galactic radio pulsars 6 X 10 7 — 4 X 10 s G. The right panel relates the 
inferred magnetic fields for the accreting systems to the accretion rate (in % of the Eddington rate). This figure indicates that the 
fields are most seriously overestimated for the fastest accreting systems. [Radio pulsar data taken from the radio pulsar catalogue 
http://www.atnf.csiro.au/research/pulsar/psrcat/ Accreting neutron star data determined from Table HI ] 



the disk. Viscosity is the main agent that dissipates energy and angular momentum, and thus enables matter to flow towards 
the central object. Since the microphysical viscosity (likely due to the mag netorotational instabi l ity in some form) is difficult 
to characterise, it is common to use the so-called a- viscosity introduced bv lShakura fc Sunvaevl Jl97ct) . i.e. let the kinematic 
viscosity be parametrised as 

c 2 

v = ac s H = a-^- (12) 

__K 

Here c s is the sound speed in the disk and H ~ c s /f_K is the vertical scale height. In this description, v is a function of r 
since both c s and £_k vary with position, but a is taken to be constant. In effect, this leads to a model where the viscosity 
ensures that the disk remains Keplerian as matter and angular momentum is transferred through the disk. 

In the case of a magnetically threaded disk, we need to provide a description of the interaction between the disk flow and 
the magnetic field. Figure [5] provides a schematic illustration of the problem. To provide a detailed model of this complicated 
physics problem is, however, not a simple task. Nevertheless, one may hope that a somewhat simplistic description will be 
able to capture the main features of the complete problem. 

From the (^-component of the Euler equations for the disk flow we can estimate the radius at which magnetic stresses 
balance the material stresses. We thus find 

M^- \^{r)r 2 ] = -r 2 B v B z (13) 
ar L J 

where the mass transfer rate M will be assumed constant throughout the disk. This relation illustrates the difficulty involved 
in constructing a consistent model. If we consider a thin accretion disk, then the z-component of the magnetic field can be 
taken to be that associated with a rotating dipole (with dipole moment /i) 

B> = -£ = -B°(f) 3 (") 

where Bo is the surface field of the star. Even though the field may be much more complicated close to the stellar surface 
the dipole contribution will dominate far away. The problem is associated with B v . This component, which vanishes in the 
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absence of a disk, represents the degree to which the magnetic field is dragged along with the matter flow. It is this interaction 
which leads to the torque on the star that we are ai ming to model . 
From the MHD induction equation we find that llMestellbooi) 

d t B v w ^ « V x (v x B) = 7 (fi - n K )B z (15) 

where the star (and the magnetic field) is rotating at the constant rate fl. I n this equat ion, it is assumed that the disk 
flow changes from quasi-rigid to Keplerian over a lengthscale -R/7, with 7 > 1. IWanj lll995F) has considered several different 
mechanisms for the timescale t v (and by implication the toroidal component of the magnetic field). He concludes that the 
various models lead to quite similar predictions for the accretion torque. This is fortunate, since it means that the model is 
not very sensitive to the unknown physics. Here we will assume that the m ain mechani sm that prevents the field from being 
dragged along with the flowing matter is turbulent diffusion. This leads to IjWanel l 1995*1. 

H 1 

v ac s cxQk 
and consequently 
7 fi - K 



(16) 



(18) 



o B * ( 17 ) 

We can now return to Eq. 1131 and determine the "inner" edge of the accretion disk Rq, at which the matter flow departs 
significantly from a Keplerian profile; 

/Ro\ 7/2 = 2N C \ 1 _(Ro\ 3/2 ~ 
\Rc) MyJGM R c [ \Rc) 

where R c is the co-rotation radius and we have defined 

m - 7 ^ - 7 r 2 r 3 - 7 r 2jR6 rial 
N "-aRj-a B ^-a B °m (19) 

Where B z is the magnetic field at R c and Bo is the field at the surface of the star as before. 

We can also account for the torque due to the magnetically threaded disk outside _Ro- The corresponding torque follows 

(essentially) from integrating Eq. (1131 and we get 
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As discussed previously, the region Rq < r < R c contributes a (positive) spin-up torque, while the region R c < r < oo provides 
a (negative) spin-down torque, cf. Figure |Sj 

Finally, assuming that the matter gives up all its angular momentum (relative to the frame of the star) upon reaching 
Ro, i.e. that the matter flows along the field lines like "beads on a wire" in the region where the magnetic field dominates the 
flow 2 , we find that the total accretion torque is 



N = MVGMRo 



1 - 



(-) 



3/2' 



+ AT disk = -MVGMRo 



7/2 - 7(R /R C ) 3/2 + 3(R /R C ) 



1 - (R /R c ) 3 / 2 



This result shows that the system reaches spin-equilibrium (N = 0) when 
'Ro\ 3/2 7-V7 



Ro « 0.8R C 



(21) 



(22) 



McJ 6 

This should be compared to the result of lWand <ll995F) . The difference arises from the fact that Wang uses Eq. J7J rather than 
Eq. for the material torque at the inner edge of the disk (now at Ro instead of Rm). 

Having added the spin-down torque exerted on the star by the outer parts of the disk we find that the system reaches 
equilibrium slightly before Ro reaches R c . Nevertheless, the predicted spin-period at equilibrium 
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(23) 



does not differ much from the more naive prediction provided by Eq. Illfl . Of course, the actual spin period at equilibrium 
now depends explicitly on the ratio a/7. Unfortunately, both these parameters are largely unknown. In addition, there are 
many uncertainties (at the level of factors of order unity) in the model. 

I n order to proceed we note that the viscosity par ameter a is usually assumed to lie in the range 0.01 — 0.3 jFrank et alJ 
2002), while 7 has been assumed to be of order unity iWandll987l) . If we consider values in this range, what does the model 
imply for the magnetic fields of the accreting LMXB neutron stars? From Eq. 1231 we find that a canonical neutron star will 
have equilibrium of 3 ms if 

V2 / »> \ 1/2 



B « 9.4 x 10 s 



M 
M Ed < 



(24) 



We see that for a/7 ~ 0.1 a star accreting at the Eddington rate is predicted to have a magnetic field within the range deduced 
for the millisecond radiopulsars. On the other hand, a star accreting at 1% of this rate will require a larger value of order 
a/7 ~ 1 in order to lie in the range indicated in Figure This means that the inclusion of the torques from a magnetically 
threaded thin disk is, in principle, sufficient to remove the direct need for an addition spin-down mechanism like gravitational 
radiation in these systems. Of course, this is achieved at the cost of introducing the poorly constrained parameters a and 7. 
If we want to adjust these parameters in such a way that the inferred magnetic fields agree with those for the radio pulsars in 
Figure we essentially need to introduce a suitable B^ = B V (M). Despite this possibility, we do not think that the thin-disk 
model is entirely satisfactory. As we will argue in the next section, additional physics should be included in order to describe 
the fastest accreting systems. In essence, this means that we will only rely on the thin disk model for systems accreting below 
a few percent of the Eddington rate. From the above estimates we see that these systems are adequately described if we take 
0/7 w 1. Hence, this will be our canonical value from now on. 



4 THICK DISKS NEAR EDDINGTON ACCRETION 

The thin disk model we have discussed so far is able to explain many features of accreting neutron star systems. Yet we will 
see that it cannot be relied upon for rapidly spinning stars accreting near the Eddington limit. Given this, it is meaningful 
to ask what the crucial missing piece of physics is. At this point, the most naive assumption in our discussion concerns the 
accretion torque arising from the inner edge of the disk, at Ro. While it seems reasonable to assume that the matter moves 
along the magnetic field lines in the inner region for low rates of accretion, it is not so clear that this model will work for 
faster accretors. Several mechanisms may alter the picture. Obvious possibilities are: radiation pressure from the emerging 
X-rays, the near balance between centrifugal and gravitational forces for rapidly spinning stars, heating of the disk in the 
inner region etcetera. 

2 In reality the problem is expected to be significantly more complicated, with the answer depending on the detailed physics in an 
extended transition region. 
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Figure 3. The main lengthscales in the accretion problem. The relevant parameters are taken to be a/7 = 1, Bo = 10 s G and the star is 
assumed spin at a rate corresponding to equilibrium for thin disk accretion. The standard magnetosphere radius (for spherical accretion) 
Rm is shown as a thin dash-dotted line, and the corresponding radius for a thin magnetically threaded disk Rq is a thin solid line. The 
co-rotation radius R c , at which a Keplerian disk co-rotates with the star, is a thick dashed line. Finally, the distance at which radiation 
pressure balances gas pressure R m i is shown as a thick solid line. The figure shows clearly that R m i >> R c > Ro > Rm above a few 
percent of the Eddington accretion rate. This suggests that radiation pressure must be accounted for, likely leading to a thickening of 
the disk and a sub-Keplerian flow in the inner region. 



As a first stab at including these effects we will consider the radiation pressure. One can show that radiation pressure 
balances the gas pressure at a radius jFrank et al ] |200i IPadmanabhanlEoOlh 

( • \ 16/21 / v 1/3 

M EddJ \ M ®J 



where 



f 



R\ 1/2 



(?) 



-, 1/4 

(26) 



Strictly speaking, this result holds only for non-magn etic disks, but one can argue that it remains a good approximation also 
in the magnetic case JCampbell fc Heptinstanl ll99&l. Moreover, it is easy to show that the factor involving / will be near 
unity apart from in the absolute vicinity of the stellar surface. Hence, we can use 



/;„„ * sWa 2/21 [ ^— } I —fij- ) km (27, 



16/21 / M X 1/3 



(^Edd) {lAM t 



as a good approximation. Let us contrast this to the standard radius of the magnetosphere, Rm- We find that R mi = Rn 
when 



r- 2xl ° a 77^ (28) 

The key lengthscales in the problem are illustrated in Fig. 01 This figure shows that, for a neutron star with a weak magnetic 
field (a typical millisecond pulsar) radiation pressure will be important for accretion rates above a few percent of the Eddington 
rate. 

We thus have to ask how the radiation pressure affects the model outlined in the previous section. Phenomenologically, 
the disk is likely to expand leading to the flow becoming sub-Keplerian. In fact, that the thin disk model is unstable in a 
region where t he radiation pressure dominates the gas pressure was demonstrated a long time ago bv lLightman fc Eardlevl 
l)l974ll (see also lshaniro et all jl97fij) V In order to account for this quantitatively, let us consider the following model. The thin 
disk description is relevant outside R m i, and hence describes systems accreting below the critical rate. For faster accretors , 
there will exist an inner region inside R m i where the disk is no longer thin. To model this region we follow lYi et alJ Jl997f) . 
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and a ssume that the flow is such that Q. = AO,k with A < 1 3 . Before moving on we should note that the study of lYi et all 
Jl997h pertains to advection dominated accretion below a critical accretion rate, while our model concerns radiation pressure 
dominated disks above a critical accretion rate. This may seem a cause for concern, especially since advection dominate d 
flows are almost exclusively used in discussions of slowly accreting systems. However, as pointed out bv lNaravan fc y1 (|l994l 'l. 
the corresponding solution to the equations describing the accretion problem is likely to be relevant also for rapid accretion. 
Furthermore, the model is sufficiently simple to serv e our present purposes. A more detailed analysis that supports the basic 
principles behind our model has been carried out bv lCampbell fc Heptinstail (Il998h . 

Repeating the arguments from the thin disk analysis, we find a new co-rotation radius A 2 ^R C and the inner edge of the 
thick disk region R' is now determined from 



R'o 



A 2 /'iR c 



7/2 



2N' 



A 4 / 3 MVGMR C 



, _ _ I 5 



A V R, 



3/2 



where we have defined 

,2 



c aA 2 R% 



= ^-B 2 Z A 2 R\ 



The torque from the inner disk region follows from 



iVthic 



M 2 7 



-7 



3/2' 



dr 



while the outer (thin) disk contributes a torque 

iVthin = - 



B^B z r 2 dr = ^ 



[°° 1 




/ r \ 3 / 2 " 




i - 


\R C ) 



dr 



Working out the algebra, we find that the total torque can be written 



N = My/GMR'o — 
where 



7 7Q 2 ACL -A) 



( Rmi ) 



3/2 



LO = 



3/2 



In this slightly more complicated model, the system will reach spin-equilibrium when 
7uj , _ 2 , A{l-A) ( R c ^ 3/2 



(— V 

V Rmi ) 







(29) 



(30) 



(31) 



(32) 



(33) 



(34) 



(35) 



6 3 3 

Since we must have R' < R m i we are always interested in the smallest of the two roots to this quadratic. The problem 
simplifies considerably if we note that 

-16/21 



Rm 



2 x W~ 2 a~ 2/21 



M 



U ms/ 



2/3 



(36) 



for a canonical neutron star, cf. Fig. [3] This means that, for a sizeable fraction of the Eddington accretion rate and millisecond 
spin periods, we have equilibrium when 

V7 



7- 



6 



R'o » 0.8A 2/3 R C 



From this we can infer the spin period at equilibrium; 



-10/7 



-3/7 



M 



-3/7 



Bp 
VlO 8 G 



6/7 



M 



I AM, 



5/7 



R 



10 km 



18/7 



(37) 



(38) 



This result differs from the thin-disk model only by the factor of A. However, it is easy to see that this is a key factor which 
may lead to considerable differences in the predicted equilibrium spin periods. 

To complete the thick disk model, we need to estimate the coefficient A which describes the nature of the sub-Kepler ian 
flow. To do this, we consider the radial Euler equation which (for a thin disk) can be approximated by jFrank et alJl2002h 



3 This model should provide an acceptable representation of the inner disk flow, but there are important caveats: All radiation dominated 
config urations tend to be subject to thermal and convective instabilities and hence may not be stationary, see eg. ISzuszkiewicz fc Miller! 
12001) for discussion. 



10 N. Andersson et al 



dv r v% id f B 2 \ GM B% 

V r — J Pl F — + ~r^ ( 39 ) 

or r p or V 87r I r A VKpr 

This equation will rem ain approximately relevant in the case of a thick disk provided that it is interpreted as a height average 
jNaravan fc Yj|l995allbV Apart from very near the Eddington accretion rate the dominant velocity component is v v . The 
situation near MEdd is complicated by the fact that the matter in the disk becomes highly virialised. In our thick disk model, 
we expect the radiation pressure to dominate in the region R' < r < R m i- (It is worth noting the difference between the radial 
and azimuthal Euler equations here. In the latter the axisymmetric radiation pressure will not play a role and the magnetic 
and viscous stresses dominate.) 

We express the radiation pressure gradient in terms of the co-moving radiation flux L co iMillerlll990l lMitralll99 Al 

rfffrad _ HP L co 

dr c 47rr 2 

where k is the opacity of the matter. Since the Eddington luminosity follows from 

A-kGMc .... 
^Edd = (41) 

we have 

dp rad GM L co 

-*r = - p —L^ (42) 

Using this relation in Eq. 13911 we see that the velocity profile becomes sub-Keplerian with 



(40) 



— where A = Jl--^- (43) 
r V ^Edd 

As a rough approximation we can assume that the co-moving flux is equal to the stationary flux observed at infinity Lx = 
GMM/r where r is the distance to the source. Then 

L °° - A (44) 



L E dd M Edd 
and we see that 



/ M 

v v w Arfl K with A = \ / 1 : (45) 

V M Edd 

(effects due to eg. beaming have obviously not been included in this estimate). 

The results we obtain by combining this approximation with the predicted equilibrium period for the thick disk model are 
illustrated in Figure |1] This figure shows that the thick disk model leads to significantly longer equilibrium spins for rapidly 
accreting systems. Conversely, we can use Eq. (|38^ to deduce a system accreting at 90% of the Eddington rate and which is 
observed to spin with a 3 ms period, should have a magnetic field of B w 1.4 x 10 8 G. A field of this strength would put 
this system well within the range of fields inferred for the millisecond radio pulsars, cf. Figure 2] The figure shows that our 
thick disk model leads to predicted magnetic fields for the LMXBs which accord well with those of the galactic millisecond 
radio pulsars. (In order to infer the magnetic fields shown in Figure 0] we have assumed that the fastest accreting systems 
have M — 0.95A/Edd- This is somewhat ad hoc, but it should be noted that the model breaks down, in the sense that A — > 
which leads to P cq diverging, as M — > M Edd . There is also significant uncertainty in the accretion rates given in Tabled) 

The fact that radiation pressure will affect acc retion disk struc t ure, a n d hence the spin period of ne utron star s in LM XBs, 
has pr e viously been discussed by seve ral authors Jwhite fc Stellal Jl988h : lGhosh fc Lambl Jl99ll . fl992^> . see also iMiller et al] 
( 1998); IPsaltis fc Chakrabartvl (1999)). These models, like ours, give lower inferred magnetic fields for high accretion rate 
sources when radiation pressure is taken into account. One issue associated with the previous models is that if one assumes 
sp in equilibrium, the models predict a strong correlation between magnetic field and accretion rate (see for example Eq. (25) 
of IMiller et all (1998')). No direct measurements of LMXB magnetic fields have yet been made, so this correlation cannot be 
tested, but the physical basis for s uch a strong relation is at best unclear. In fact, this has been one of the arguments against 
magnetic spin equilibrium models llBildstenlll998l : Tushomirskv et alJl2000l) . The model outlined in this paper also predicts a 
correlation between magnetic field and accretion rate. The "dependence" of B on M is however weaker, due to the dependence 
on accretion rate of the factor A. This illustrates that small modifications to the accretion model may be able to remove some 
of the perceived difficulties associated with magnetic equilibrium models. 

We conclude this section with a brief discussion of the observational consequences of this model with regard to the 
detection of X-ra y pulsars. Naively one expects the X-ray pulsars to have higher inferred magnetic fields than the other, non- 
pulsing, LMXBs iCumming et al.ll200ll) . As is clear from Figure|l| the thick disk model does not lead to the pulsars clustering 
at higher magnetic fields than the other sources. One possibility, suggested bv lTitarchuk et all (120021) . is that we are prevented 
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Figure 4. The predicted spin periods at equilibrium for the thick disk model, for a = ■y = 1. In the left panel we show P cq as function of 
the accretion rate for magnetic fields which bracket the range for the millisecond radio pulsars: B = 6 X 10 7 G (thick dashed curve) and 
_B = 4 X 10 s G (thick solid curve) . For comparison we also show the prediction of the naive model where spinup ceases at R c = Rm (thin 
solid curve). The right panel compares the inferred magnetic fields for LMXBs to those of the radio pulsars and should be compared to 
the right panel in Figure 1. 



from seeing pulsations in many systems due to atmospheric scattering. A preliminary study bv lKrauss fc Chakrabartvl \2QQ- 
suggests that the scattering hypothesis may not be borne out by the data, but this is an area of ongoing research. 



5 CONCLUSIONS 

We have discussed the accretion spin-equilibrium for neutron stars in LMXBs. The outcome of this study is a more detailed 
model of the accretion torques and an appreciation that it is possible to construct a reasonably simple and consistent model 
for these systems without invoking additional spin-down torques due to, for ex ample, gravitationa l radiation. This result is not 
particularly surprising. After all, the accretion problem is extremely complex iFrank e t alj|20 02t). and the torq ues c onsidered 
in the studies that argued for the need for an additional spin-down mechanism (see Iwhite"^; Zhand (Il997l) and iBildster] 
(1998)) were somewhat simplistic. 

Of course, our results should not be taken as proof that the LMXBs do not radiate gravitational waves. The various 
proposed mechanisms for generating asymmetries in rapidly spinning, accreting neutron stars remain (essentially) as viable 
as before. The key difference is that we have eliminated the rationale for locking the gravitational radiation luminosity to 
the non-magnetic torque Ms/ GMR, which has been used as an order of magnitude estimate in most studies to date. In our 
picture, one would not be able to infer how the spin down due to gravitational radiation combines with the accretion torque 
from the observed spin periods. This alleviates s ome "problems" with the gravitational-wave models. In the case of accretion 
induced asymmetries in the crust feildstenll 998). one can show that the quadrupole deformation required to balance accretion 
is 

1,1 f — ) 1 (^y 5/2 



V^Edd/ V3ms, 

This s hould be compared to the maximum deformation that the crust can sustain, which according to llJshomirskv et, ail 
J2OO0I) can be approximated as 



< 5 x 10 



-7 / ttbroak \ 



(47) 



where the breaking strain itbrcak is usually (based on results for terrestrial materials) assumed to be in the range 

These estimates show that the breaking strain must be near the upper limit of the expected range in order for these asymmetries 

to balance near Eddington accretion in a star spinning at a period of a few milliseconds. By weakening the accretion torque, 

while at the same time not altering the mechanism generating the asymmetry (eg. the accretion rate), this issue is made less 

critical. 
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Our results also impacts the suggestion that the gravitational waves are emitted by unstable r-mode oscillations in the 
stellar fluid. In this case, the r-modes are expected to become unstable below a critical rotation period P cr it. The point at 
which the i nstab ility becomes relevant depends on many complicated issues concerning viscosity, superfluidity etcetera (see 
lAnderssonl (|2003h for a discussion) but it is plausible that P cr it ~ 2 — 3 ms. In the context of the present model, we obviously 
need Pcrit > P eq in order for the r-mode instability to be relevant. Considering the results illustrated in the left panel of 
Figure [4] we expect that the instability may come into operation in weak magnetic field systems which are neither very slow 
nor very fast accretors. 

The most important next step in modelling the LMXBs con cerns the variability in the spin with varying accretion rate. 
The spin of accreting X-ray pulsars is known to vary considerably JlBildsten et alll997h on a timescale which is roughly similar 
to the variations in the accretion rate. But the data also suggests that there may not be a direct link between increased X- 
ray flux and an increase in the spin-up torque. It is important to understand this variability in general. This is also a very 
important issue for attempts to search for gravitational waves from the LMXBs. Any variability on timescales shorter than 
the observation time that remains unaccounted for will likely lead to a significant loss in signal-to-noise ratio. Our aim is to 
turn our attention to this challenging problem in the near future. 
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